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Abstract 

An infinite log-gas formalism, due to Dyson, and independently Fogler and Shklovskii, 
\ is applied to the computation of conditioned gap probabilities at the hard and soft 

^ [ edges of random matrix /3-ensembles. The conditioning is that there are n eigenvalues 

in the gap, with n <^ \t\, t denoting the end point of the gap. It is found that the 
entropy term in the formalism must be replaced by a term involving the potential 
drop to obtain results consist with known asymptotic expansions in the case n = 0. 
With this modification made for general n, the derived expansions — which are for the 
logarithm of the gap probabilities — are conjectured to be correct up to and including 
terms 0(log|i|). They are shown to satisfy various consistency conditions, including 
an asymptotic duality formula relating /3 to 4//3. 

1 Introduction 

Fundamental to Dyson's [13j pioneering studies of random matrices is the log-gas picture. 
The essential point is that the joint eigenvalue probability density function for various 
random matrix ensembles can be computed exactly, and is seen to be identical to the 
Boltzmann factor for certain classical particle systems interacting via a repulsive logarith- 
mic pair potential. These so-called log-gas systems are in equilibrium at particular values 
of the dimensionless inverse temperature /?, restricted to the value /? = 1 in the case of the 
matrix ensemble exhibiting an orthogonal symmetry, (3 = 2 for a unitary symmetry and 
/3 = 4 for a unitary symplectic symmetry. In the application of random matrix ensembles 
to quantum mechanics these cases correspond to a time reversal symmetry T with = 1, 
no time reversal symmetry, and a time reversal symmetry T with = —1, respectively. 

The significance of the log-gas picture is that it permits the use of macroscopic argu- 
ments to compute certain probabilistic quantities in the bulk scaling limit (where the num- 
ber of eigenvalues goes to infinity with their bulk density fixed to the value p). Specifically, 
the particle system is approximated as a continuum fluid conflned to a one-dimensional 
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domain but obeying the laws of two-dimensional electrostatics. As an example, consider 
£;^""'(0; (0,t)), the probability that there are no eigenvalues in the interval (0, t) for the 
bulk scaling limit of a random matrix ensemble indexed by /3. Dyson's first application of 
the log-gas analogy was to derive the asymptotic form 

log^,^"^>^(0; (0,t))^^^ -/3i^ + (f - ? + ^'-'^ 

Although it was remarked [13\ above eq. (92)] that the variational principle underlying 
(jl.ip may contain errors of order log pt, other arguments were presented that suggested the 
log-gas prediction of the explicit term 

^i^.o.P«. (1.2) 

may in fact be correct at this order. 

A noteworthy feature is that (jl.ip is valid for all /3 > 0, and thus in particular applies 
uniformly to all three random matrix symmetries. In fact, the general /3 > form of 
(|l.ip is relevant to contemporary random matrix theory [23j. Thus there are random 
tridiagonal matrices jlT] and random Hessenberg unitary matrices [23] that /3-generalize 
the eigenvalue probability density functions for the Gaussian and circular ensembles from 
classical random matrix theory. A remarkable development relating to these constructions 
has been the characterization of £'^"^'^(0; (0, t)) in terms of a stochastic differential equation 
[25] . [33] , This was subsequently used [3l] to give a rigorous derivation of (jl.ip . and to 
furthermore establish that the explicit form of the O(logpt) term is not (II. 2[) but rather 

^(f + |-3)logpt. (1.3) 

Generalizing £^^'^^'^(0; (0,t)) is the probability £'^""^(n; (0, t)) of there being exactly n 
eigenvalues in the interval {0,t) for the bulk scaling limit of a random matrix ensemble 
indexed by f3. Many years after his pioneering work f 13] . Dyson [H] and independently 
Fogler and Shklovskii [TS] used the log-gas picture — now applied directly to the infinite 
bulk state (in [13] use was made of a scaling limit of a large deviation formula in the finite 
system; see ^ — to derive the generalization of (jl.ip 



+ Sfl-f-^Vlog^ + l). (1.4) 



2V 2 2 J \ n 



In fact, in the cases /3 = 1, 2 and 4 this same formula, except requiring n to be fixed, had 
been obtained earlier [2\ using an analysis based on explicit Painleve/Fredholm determinant 
expression for E^'^^^{0; (0,t)). 
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The primary aim of the present paper is to extend Dyson's p3j log-gas study of the bulk 
conditioned gap probability £'^""^(n; (0,t)) to the analogous quantities £'^°^*(n; (t, oo)) and 
^hard^j^. (Q, i)) Corresponding to there being exactly n eigenvalues in the intervals (t,oo) 
and (0,t) about the soft and hard spectrum edges respectively, for the regime 1 <C n <C 
|t| — ?■ oo. The soft edge refers to the neighbourhood of the largest or smallest eigenvalue, in 
the situation that the corresponding eigenvalue density exhibits a square root profile, which 
from our formalism extending [H] is taken to be (semi) infinite in extent. The hard edge is 
the neighbourhood of the smallest eigenvalue for matrices with non-negative eigenvalues, 
when the spectral density exhibits an inverse square root singularity profile, which again 
is to be taken of (semi) infinite extent. 

As examples of explicit realizations of these edge states in terms of limits of finite N 
matrix ensembles, let ME^^Ar(g') denote the eigenvalue probability density proportional to 

N 

n^cAo n iAfc-A,i^. 

1=1 l<j<k<N 

The choice g{X) = e~^^ 1'^ defines the Gaussian /3-ensemble, and the choice (7(A) = 
^/3a/2g-/3A/2 > q) defines the La guerre /3-ensemble. We know from [19j that the largest 
eigenvalues are, respectively, to leading order in at A = \/2N and A = 4A^, and that 
[TUl ES] upon the scalings 



Az = ^/2]V + -^, A.=4iV + 2(2iV)V2x, 

the correlations have the same well defined — )• 00. The corresponding scaled density 
admits an explicit /3-dimensional integral formula in the case of /3 even [10], and 
this exhibits the square root singularity characterizing the soft edge state by way of the 
asymptotic formula 



C''(x) ~ (1.5) 



The standard realization of the hard edge is the Laguerre /3-ensemble, upon the scaling 

[EKE] 

\ - 



For a = c — 2//3, ca positive integer, and /3 even the corresponding scaled density p^^^^''^ 
admits an explicit /3-dimensional integral formula [20j . Asymptotic analysis of this integral 
shows 

pf^.^'^ix) ~ TT^, (1-6) 



X 
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thus exhibiting the characterizing density profile at the hard edge. 

The relevance of the spectrum edge to problems in theoretical and mathematical physics 
was not apparent in the era of Dyson's first works. In fact, at the soft edge there is a 
somewhat hidden, but nonetheless fundamental, relationship to growth problems of the 
KPZ class (see e.g. the recent review [IZ]). And in the case of the hard edge the statistical 
properties of the smallest eigenvalue can be related to data from lattice QCD simulations 



A priori the log-gas formalism can be applied to the study of the asymptotic form of 
gap probabilities in the cases the gap contains no eigenvalues, or in the cases there are n 
eigenvalues with 1 <C n <C In the former case the asymptotic form of the logarithm 
of the gap probability for general (3 at both the hard and soft edges is known from other 
considerations |20t [5l HI [29], up to and including terms 0(log|i|). This information is 
important for two reasons. First, it tells us the infinite log-gas formalism at both the hard 
and soft edge is in error in the second order (free energy/ entropy) contribution. But 
fortunately, guided by the bulk case, an ansatz can be made to correct this error. Second, 
it allows the asymptotic expansions obtained for 1 <C n <C |t| to be very naturally extended 
to hold (as a conjecture) uniformly in n, up to and including terms 0(log \ t\). 

These latter expansions, which are the main result of the paper, read 



As a check, using results from [22], asymptotic duality equations — relating gap probabil- 
ities with given /3, n to gap probabilities at 4//3 and with modified n and length scale — 
are identified for both S|^^^'^(n; (0, t)) and Ef^^{n]{t,oo)), and pTT]) and (fL8]) are shown 
to exhibit the required inter-relations. 

2 Strategy of the log-gas argument 

2.1 Formalism for a general background potential 

The finite log-gas system consists of N mobile charges, of strength +1 say, and a smeared 
out neutralizing background charge density —pf,(x). We are interested in the N ^ oo limit. 



(see e.g. [32]). 




and 




log\t\-^/\ (1.8) 
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in the circumstance that Pb{x) is supported on some semi- infinite or infinite domain /. A 
physical argument tehs us that an equihbrium state will exhibit local charge neutrality, 
and so = Pb{x), where denotes the partial density (otherwise the charge 

imbalance would create a non-zero electric field, which in turn would mean the system is 
not in equilibrium). 

Let Ej3{n; J)[pb] denote the probability that in this system there are precisely n eigen- 
values in the interval J. The basic hypothesis (see e.g. [231 §14.6.2]) underlying the log-gas 
argument is that for | J| — )■ oo 

logEp{n;J)[pb] = -PSF, (2.1) 

where 6F is the change in energy resulting from conditioning the equilibrium particle 
density so that the interval J contains n particles. Furthermore, the change in energy 
consists of two distinct parts, namely the electrostatic energy Vi and the free energy V2, 
the latter being entirely due to entropy. 

The electromagnetic energy is responsible for the leading order (in | J|) contribution to 
(12. ip . With /9(i)(x) now denoting the new equilibrium particle density, and 



(Pix) = - J^ip(i)iy) - Pbiy))^og\x - y\dy, (2.2) 

denoting the electrostatic potential at a point x due to the perturbed charge density, we 
have 

Vi = ^ j{p(i) (x) - pb{x))4>{x)dx. (2.3) 

The second-order contribution to (j2.ip results from both the electrostatic energy and the 
free energy. Let fp = //3[/0(i)(x)] denote the density function for the free energy per particle 
at point X. The change in this density function for the free energy per particle due to the 
perturbation is 

1 1\, ( P(i){x) 



fp[P{i){x)] - fp[pb{x)] = Q - ^) log 



Pb{x) 



where the equality follows from the exact expression for //3[/o] (see e.g. [23] §4.8.1]). Thus 
the total change in free energy is 

With the background charge density Pb{x) given, two-dimensional electrostatics is to 
be used to determine the complex electric field 

Eiz) = f 'Jn^^l^^d., (2.5) 

It Z-X 
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which in turn is used to compute the particle density and electrostatic potential 

4>{x). The validity of this approach is restricted to large gap size \J\. Thus it is only in this 
circumstance that we are able to use physical principles to fix (p{x) in certain domains and, 
furthermore, to specify some properties of to the extent that E{z) can be uniquely 

determined. 

To see this suppose for definiteness that 1 <^ n <^ \ J\. We then expect the n particles 
to form a macroscopic blob within some interval J' C J, and thus 

/9(i)(x)=0, forxGJ\J', (2.6) 

while 

P(i)(x)(ix = n. (2-7) 

J' 

Furthermore, we expect (x) to vanish as x approaches the boundary of J' from within J' 
(a soft edge singularity), while we expect /0(i)(2;) to diverge as x approaches the boundary 
of J from within J (a hard edge singularity). It could also happen that the left hand end 
points of J and J' coincide, in which case a hard edge singularity of /9(i) (x) is expected. 

From the viewpoint of macroscopic electrostatics the regions J' and I\J, being occupied 
by what is effectively a continuum of mobile charges, behave as perfect conductors and thus 
have (pix) equal to a constant throughout. As a normalization, in the (semi) infinite region 
I\J this constant can be chosen to equal zero. Then for x E J' we will have (pix) = —v, 
with V > 0. Furthermore, (p{x) must vary continuously from this value to zero, and, with 
t the smallest end point of I\J and the largest end point of J', 

^dx = V. (2.8) 
T' ax 

The fact that the particle density vanishes in the complement of the regions of constant 
potential allows (12. 3p to be simplified to read 

^1 = ^ / {p(i){x) - Pb{x))dx I pb{x)(j){x)dx 
^ J J' ^ Jj\J' 

= -'^-\j^Pb{x)mdx, (2.9) 

where to obtain the second equality use has been made of (j2.7p . And with the left hand 
end point of J denoted Jq, integration by parts together with the fact that (l){Ja) = can 
be used to conclude from ()2.9p that 

v. = ~AI ( rMy)<iy)fi^. (2.10) 



2 2 Jj\j, \Jj^ J dx 



Since pb{y) is given, we see from (I2.10p and ()2.8p that Vi depends only on and not 
on both P(i)(a;) and 0(x) as appears in (j2.3p . But knowledge of Pi^i){x) is required for the 
computation of V2 as specified by (12.41) . 
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2.2 The background density pb{x) = ex" with n = 

We now specialize to a semi-infinite domain / = (0,oo) and background density Pb{x) = 
cx'^,a > —1. This includes as special cases the soft edge (a = 1/2; recall (jl.Sp ) and 
hard edge (a = —1/2; recall (|1.6p ). It also includes the density profiles of the so-called 
multicritical matrix ensembles, for which a = 2k + 1/2, k G Z>o; see e.g. p)j. Let J = (0, t) 
be conditioned to have precisely n eigenvalues. Although in the above discussion it has 
been assumed n » 1, repeating the reasoning for the case n = shows that the appropriate 
limit of (j2.10p remains valid, and thus 



which is, of course, a property of pb{x). Assuming (I2.12p . the complex electric field ()2.5p 
exhibits an analogous scaling 




(2.11) 



In this special case the particle density p(^)(x) depends only on x and t. Moreover, as 
observed in [19], the t dependence can be factorized by the scaling ansatz 




(2.12) 




But 



dx 



lim ReE{z) 



(2.13) 



and so 



dx 




dx 



t=i 



(2.14) 



Making use of this latter formula in (|2.1ip we obtain 




(2.15) 



and thus, upon substituting in (|2.ip for 5F, the prediction [TS] that for t — oo 



log£;^(n;(0,t))[cx°] ~ -/3ct2("+i). 



(2.16) 
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To go beyond this scaling prediction, and compute for example the explicit form of c 
in (|2.15|) . we seek to deduce the explicit functional form of E{z) as a consequence of its 
necessary properties. In relation to the latter, we begin by noting that the charge neutrality 
condition 

j{p[i){^) - Pb{y))dy = 0, 

implies 

E{z) ~ o(z-i). (2.17) 
Another condition is that as t — )• oo we require 



Eiz) -c / dx, (2.18) 

Jo z-x 

which is the electric field due to the background (this is well defined for — 1 < a < 0; for 
other a > it is specified by analytic continuation). 
By considering the contour integral 

aw. 



Y.Z-W 



where S is the contour traversing the ray 9 = argz from infinity down to the origin, 
indented at the point z, then the positive real axis, we can evaluate (|2.18p to obtain 



E(z) -cz'^ PV / dx + vricz" (2.19) 

1 — X 







= (^-PV j - — -dx + TTiJ , (2.20) 

where PV denotes the principal value. By breaking the integral at x = 1 and changing 
variables x ^-^ 1/x in the range x G (1, oo) the principal value integral can be evaluated to 
give 

^(z) ^ cz"(-7rcot7ra + 7ri). (2.21) 

We see immediately that the values a = {k — 1/2), /c = 0, 1, . . . are distinguished, as then 
E{z) is pure imaginary, and 

E{z) TTicz°' as t oo. (2.22) 

This turns out to be crucial in constructing E{z) so that, for example, (j2.13p vanishes for 
X > t, and so restriction to these values of a will be assumed henceforth. 
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To continue, we use the fact that the only singularity of E(z), apart from z = as 
already present in Pbiz), is z = t. Furthermore, this is required to be an inverse square 
root, corresponding to an (induced) hard wall singularity. We see immediately that for 
Q = -1/2 



E(.) = „c-L (l - yCII) , (2.23) 

has this property, as well as the necessary properties (|2.22p and (|2.17p . Furthermore, the 
real part of E{z) on the real axis vanishes for x > t, telling us via (|2.13|) that the potential 
is a constant in this region. And the general formula 

P(i){x) - Pbix) = -- lim Im E{z), (2.24) 

which follows from ()2.5p when applied to p.23p gives p(i)(x) = for Q < x <t. Hence, we 
conclude that (j2.23p satisfies all necessary conditions required from physical principles. 

We now proceed to use (I2.23P to compute Vi. Thus substituting in (12.13P gives that 
for < X < t 

d(b vrc 



dx ^t- X 

Substituting in (j2.1ip , scaling the t dependence and computing the integral (a special case 
of the Euler beta integral) gives 

Vx = (2.25) 

thus specifying c in the case a = —1/2. 

For a = k — 1/2 with k G we generalize ()2.23p to read 



E{z) = vricz" ( 1 - J-—pk{tz-') ] , (2.26) 



where 



j=0 



x=0 



(^),:=n(^ + l)...(n + j-l). 

3=0 ^' 
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This modification is needed for the necessary condition (j2.17p to be satisfied, while leaving 
unchanged the validity of the other necessary conditions. It tells us that for < x < t 



^ = Trc^£=pk{tz-^), (2.27) 
ax — t 



and this in turn substituted in (|2.1ip implies 

^ ^[£L,2fc+i V r(2fc-j + 3/2)r(l/2) (-1/2), 
' 2A; + 1 r(2A;-j + 2) j! ' ^ " ^ 

The simplest case is fc = 1, corresponding to ph{x) = C\/x, for which ()2.28p simplifies to 
read 

= (2.29) 

In the cases a = ±1/2, the background density Pfe(x) = ex" has distinguished values of 
the proportionality constant 

c=|f' ""7^. (2.30) 

These are read off from (|1.6p and (jl.Sp respectively. Substituting the values of c from (j2.30p 
in (|2.25p and (j2.29p respectively, then substituting the corresponding values of Vi for 6F 



in (12. ip with n = predicts that 



iog£;f*(0;(i,oo))^^~^-^|f, (2.31) 



log4-'i(0;(0,t)) ~ -'-t. (2.32) 



/5 



The first of these predictions agrees with results of recent rigorous asymptotic analysis 
|28] based upon the characterization of i?^°^*(0; (t, oo)) in terms of stochastic differential 



equations [16j; the second agrees with the result obtained from exact asymptotic analysis 
of an m-dimensional integral formula available f or ■E^ '^'''^(0; (0,t)) in the case that /3a/2 = 
m G Z>o,/3 > general [20]. We remark that (j2.32p has in fact been derived earlier in 
a log-gas analysis [5j. The implementation though was different to that used here, as it 
was applied to the hard edge of the finite N Laguerre /3-ensemble (see ^ below for further 
comments on finite strategies). 

In the next two sections the log-gas argument leading to (I2.3ip and ()2.32p will be refined 
to obtain the analogue of Dyson's asymptotic formulas for the bulk spacing (jl.ip and ()1.4p . 
for the soft and hard edges. We will begin with the latter, as it is then possible to make 
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use of results already obtained in [14J. But before we proceed we remark that in the case 
of exponents a = 2k + 1/2 (corresponding to the multicritical matrix ensembles) and with 
/3 = 2, the asymptotic form ()2.16p has recently been rigorously derived in 0, together with 
an explicit evaluation of the proportionality. However the value of the latter as implied by 
(j2.28p cannot readily be checked as the workings of [6] do not make explicit the value of c 
in the effective background density pb{x) = cx". 



3 The hard edge 
3.1 Direct calculation 

In the log-gas formalism the hard edge is specified by the background density 

1 



Pb{x) 



27rVx' 



and, furthermore, by an electrostatic coupling between a fixed charge of strength a' at the 
origin, and the charge density, 

POO 

V[ := -a / (p(i)(y) - Pbiy)) logy dy, (3.1) 
Jo 

where 

, a — 1 1 , , 

a' = ^ + ^, (3.2) 

(see [23]). This additional term accounts for the parameter a in the Laguerre weight 
rj.ai3/2^-l3x/2 ^ which in the hard edge scaling is independent of A^. 
Comparing ()3.1I) with (12. 2p shows 

Vl = a()){0). (3.3) 

Conditioning so that exactly n eigenvalues lie in the interval (0, t), and with 1 ^ n ^ t we 
expect the n particles to occupy J' = (0,6), 6 <^ t. This region will behave like a perfect 
conductor, and consequently 4'{x) = —v,v > and thus constant for x € (0,6). It follows 
from this that (13.31) can be written 



Vl = -a'v. (3.4) 
Furthermore, in this setting (I2.10p reads 

vn 1 ^ dS , , ^ 
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and we have too that 



-- dx. (3.6) 

dx 

Our immediate task is therefore to compute dtp/dx for x G (6, t). 

For this purpose we must generahze (12.231) to the present setting of a new soft edge 
singularity at z = b, while still upholding the other necessary conditions required of E{z), 
as discussed above (|2.23p . We see that the simple variation of (j2.23p 



has these properties. This furthermore satisfies the requirement that its real part vanishes 
on the real axis for x G (0, b) U {t, oo), as required by these regions being perfect conductors, 
while ()2.24p shows that /9(i)(x) = for x £ {b,t). Thus all physical requirements are 
satisfied, so we can now use (|3.7p for predictive purposes. 

In particular, for the task at hand, substituting (13. 7p in (j2.13p allows us to conclude 
that ior b < X < t, 



1 



dx 2 y x{t — x) 
Substituting in (j3.6p and changing variables gives 



(3.8) 



ivb y t-x 

while substituting in (|3.5p and simplifying gives 



^ - b 

2' (3-9) 



vn \ , , ^ 

Vi = -— + -{t-b). (3.10) 

The integral ()3.9p . while not an elementary function of t,b as in (j3.10p . can be expressed 
in terms of the elliptic integrals E',K' corresponding to the elliptic modulus 



k = ^Jb/t. (3.11) 

Thus we read off from [Ml eq. (2.32)] that 

v = Vt[E' -k'^K']. (3.12) 

Insight into why elliptic functions appear in so-called 'two-cut' (here (0,6) U (t,oo)) log- 
potential electrostatic problems can be found in e.g. [3] 
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The endpoint b is determined by setting J' = (0, b) in (j2.7p and thus requiring 

rb 



f 

Jo 



P(^i^{x)dx = n. (3.13) 
To determine /9(i) (x) we substitute (|3.7p in (j2.24p and read off that for < x < 6, 



, ^ 1 \ b — X , ^ 

Substituting in (j3.13p and changing variables x ^ x'^ allows, upon making use of [141 
eq. (2.31), (2.33)], the integral to be written in terms of elliptic integrals E,K of elliptic 
modulus (j3.1ip . giving 

n = — \E-{l-k^)K]. (3.15) 
vr 

From the elliptic function forms (|3.12p and (j3.15p . and with knowledge of the expansions 
in [m eq. (2.36), (2.37)], we deduce that for t oo, and with 1 < n < t^/^ and A; < 1, 

6^ = 4tV2„ _ + O(nVt) (3.16) 
v = t^/'^ -nlogt^/'^ + 0{n). (3.17) 

Here, in obtaining ()3.17p . use has been made of (|3.16p . Substituting (j3.17p in ()3.4p and 
(|3.10p then shows that in this same limit 

Vi+Vi = -- t^/\n + a') + (nV2 + a'n) logt^/^ + 0(71^). (3.18) 
8 

It remains to compute (j2.4p . In addition to the explicit form of the density in (0, 6), the 
explicit form in its other non-zero domain (t, 00) is also required. But this is immediate 
from (j3.7p and (j2.24p . which give 



and consequently we have 

(3.19) 
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(this form follows upon the change of variables x i— )• x^). An equivalent expression to this 
arises in Dyson's log-gas computation of the bulk spacing [14] . where it was evaluated in 
terms of the potential drop to give for the present expression 



V9. 



1 1\ V 



(3.20) 



Its asymptotic form follows immediately upon substituting (|3.17p . 

The simplest case is n = 0. With this assumed, adding the asymptotic form of (|3.20p 
to p.lSp and recalling (|3.ip shows 

Substituting this for 6F in (|2.ip then gives the prediction 



log4'^'-<i(0;(0,t)),~ -/3 



(3.21) 



with k = \ ~ ^ ■ But this contradicts the elementary exact result of [15p20j stating 

logi?hard(0; (0,t))U=o = -/3t/8. 

It turns out that consistency with known exact results for n = can be obtained if we 
replace ([3:20]) by 

y2=[\- ^ ) V. (3.22) 



This is precisely the evaluation of V2 obtained from its definition (|2.4p found by Dyson [14J 
in the bulk. Thus we then get (j3.2ip with k = |, in agreement with known exact results 
from [20[ [5l [29] . Hypothesizing ([3.22p for general n then gives 



Vi + Vi + V2 



t 



V^fr^+;^) + (^ + ^)logt^/2 + 0(n^ 



an\ 



2 2 J 

Substituting this for 5F in (12. ip . we then obtain the prediction that 



logE^-'^(n;(0,t)) ~ 



One check on ([3.23P is that it is consistent with the inter-relation 



t— >oo 

which in turn follows from [21^ eq. (3.5)]. 



an- >a— 1 



ai~>a— 1 



(3.24) 
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A noteworthy feature of (j3.23p is that for /3 = 2 it exhibits precisely the same functional 
form as the t — )• oo expansion of log £'2'^^'^(n; (0,t)) with n fixed [31], [231 eq. (9.89)] up to 
terms O(logt) independent of n. In the case n = the exact form of the term O(logt) 
in the t — t- oo expansion of log^^^'''^(n; (0,t)) is known for general /3 fSHl [3 [29] . Adding 
this to (|3.23p suggests that uniformly in n, for n <^ t, the expansion (jl.7p recorded in the 
Introduction is valid. 

The generalized form (|1.7p can be checked against the large t duality formula 



E^^^'^in; (0,i/s») ~ E^jf + 1) - 1; {0,t/~s,/p)) 



t—>OD 



(3.25) 

aM>/3a/2-/3+2 



where sp is an arbitrary length scale that satisfies S4/^(/3/2)'^ = s^. This asymptotic duality 
is a consequence of the exact duality 122, eq. (5.10)], proved for j3 even. The asymptotic 
form p.7p can readily be checked to be consistent with p.25p . 



We have seen in the above calculation of the hard edge gap probability E^^^^{n] (0, t)) that 



3.2 Relationship to Dyson's log-gas calculation of E^^^^{n] (0,t)) 

expressions obtained in Dyson's log-gas calculation of the bulk gap probability E^^^^{n\ (0, t)) 
appear. This can be anticipated by noting that under the change of variables A; i— )• Ap the 
Laguerre /3-ensemble ME^^7v('^''"'^^e~^'^/^) transforms to the chiral /3-ensemble specified by 
a probability density function proportional to 

1=1 l<j<k<N 

The latter has the interpretation [231 Prop 3.1.4] of a log-gas confined to the half-line 
A; > 0, but with image charges of the same sign for A; < 0, and a fixed charge of strength 
(a — l)/2 at the origin. 

In the scaled limit of this model the density at the edge x = 0, as for the bulk, 
is a constant to leading order, which here takes on the value l/vr (see ^3l eq. (7.76)]). 
Consequently, after conditioning so that there are n eigenvalues in (0,6), the task is to 
compute both (j){x) and p(i)(x) in the equation 



(j)[x) = - (p(i)(y) - l/vr) log Jx^ - y^l dy. 
Jo 

But defining /0(i)(— y) = P{i){y)-, this can be rewritten 

<l^{x) = - / (P(i) iv) - l/vr) log Jx^ - y] dy. (3.26) 
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We observe that with i— )• x, (j3.26p is precisely the equation for the potential in the 
bulk log-gas with background charge density — 1/vr, after condition so that there are 2n 
eigenvalues in (— This then explains, upon the change of variables x i— )• x^, the 
appearance of quantities from |14) in our calculation above. 



4 The soft edge 

4.1 The case /3 = 2 

In our log-gas formalism, we are specifying the soft edge by 



Phix) = , X > 0. 

TT 

In the realization of the soft edge as the scaled neighbourhood of the largest eigenvalue of 
the Gaussian and Laguerre /3-ensembles as discussed in the Introduction, this corresponds 
to shifting the origin to the location of the largest eigenvalue, and then changing the sense 
of direction by the mapping x i— )• —x. Thus x > now corresponds to the region of the 
eigenvalue support. With the system conditioned so that exactly n eigenvalues lie in the 
interval (0, i), for 1 ^ n <C t, in contrast to the hard wall case we now expect these 
eigenvalues to have support (61,62) with 61 > 0, and both endpoints 61 and 62 exhibiting 
soft edge singularities. We remark that this scenario has recently been exhibited in Monte 
Carlo simulations of the eigenvalue profile for the constrained Gaussian /^-ensemble — the 
constraint being that a fraction c > 1/2 of the eigenvalues are positive; see Fig.l for 
X < 0]. Like the region x > t, which we take to have zero potential, this region will behave 
like a perfect conductor and thus we will have (pix) = —v, v > for x G (61, 62). 
In this setting (|2.10p reads 

where 

v = [ ^dx. (4.2) 
Jb2 dx 

To compute d(p/dx we must generalize the a = 1/2, n = 1 case of (j2.26p . With 61 and 62 
now new soft edge singularities, we trial 



£(.)=i(yi-:^i^£E5). ,4.3) 

Observe that for this to fall off faster than for — 00 as required by ()2.17p . we must 
have 

t = 6i + 62. (4.4) 
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Furthermore, the necessary property (j2.22|) (with c = l/7r,a = 1/2) is satisfied, and z = t 
is a hard edge singularity as required. The real part of (|4.3p vanishes on the real axis for 
X = (61,62) U {t,oo), in keeping with these regions being perfect conductors, and use of 
()2.24p shows that /0(i)(a;) = for x G (0, 5i) U (62, ^ which again is in keeping with the set 
up of the problem. Finally, we observe that when n = 0, and thus bi = 62 = t/2 (the latter 
equality follows from (14.41) ). (14.31) reduces to (I2.26P in the case a = l/2,n = l. Thus (|4.3p 
passes all consistency checks, and so is presumed correct in general. 
We now make use of (j4.3p by substitution in (j2.13p to deduce that 



dx 



(bi - x){b2 - x) 

t — x 

{x - bi){x - 62) 
t — X 



and by substituting in (|2.24p to deduce that 



1 {x-b,){b2-x) 

P{i){x) = -1 



(0,61) 

X G {b2,t) 
61 < X < 62- 



(4.5) 



vr V t — X 

If we now substitute (j4.6p in ()2.7p with J' = (61, 62) and set 

bi = t/2 -d, 62 = t/2 + d, 

so as to satisfy ()4.4p . we see that d is determined by the requirement that 

rb2 



1 



dx 



{x - bi){b2 - x) 



Jbi V t-x 
This integral has the elliptic integral evaluation 



n. 



n 



— Jb2 \tE' - 2bxK'\ 
37r ^ L i J 



(4.6) 



(4.7) 



(4.8) 



(4.9) 



For d/t this means that ~ 1 and the small argument expansions of -E', K' give that 
in this limit 

d^ ~ \/2tn. (4.10) 
Next, substituting the second case of (|4.5p in ()4.2p 

ft 



V := ct>{t) - 4>{b2) 
This has the evaluation 



dx 



62 



{x - bi){x - 62) 
t — X 



-V62 [tE-{b2-bi)K\ 



02 



(4.11) 



(4.12) 
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We seek the expansion of this integral for fc^ — )• 1 or d — )• 0, and making use of the 
standard asymptotic formulae we have 



3 



1 + log 



16t" 



2V2t 



(4.13) 



Applying this, and making use too of ()4.10p . we obtain our sought expansion 



V ^ 



4 ft\ 



3 \2 



3/2 



, 3 V2tn , 1 
l + -^^log 



2 



t|3/4 



(4.14) 



Presently we have available the leading form of d, which according to (j4.7p parameterizes 
the endpoints of the support of the n eigenvalues, in terms of t (the endpoint of the 
prescribed interval) and n, as well as the leading two terms of the potential —v in {bi, 62); 
in terms of t and n. We next consider the electrostatic potential Vi. Substituting (j4.5p in 
(liTTl) shows 



vn t 



t/2-d 



dx {t — 2x)a 



' 2 ^ 671 



dx (1 — x)^ 



'{t/2-x)^-d^ 
x{t — x) ' 

/(I -x)2 - {2d/tf 



(4.15) 



x{2 - x) 

As in the case of the hard edge (|3.10p . the integral in (j4.15p can be evaluated explicitly. 
Lemma 1. Let 



H{u) -.= 2 dx (1 -x). 



(1 — x)2 — 
x{2 — x) 



We have 



Proof: Changing variables x 1— )• 1 — x, then changing variables = y we obtain 

-1 ' 



(4.16) 



H{u) 



i-y 



dy. 



Changing variables y ^ y + in this allows the u dependence to be scaled to give 



-dy, 
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and KT6\i follows. □ 
Using the result of Lemma [T] in (|4.15|) gives the evaluation 

vn ( 4cP 



Substituting (j4.10p and ()4.14p it follows from this that for 1 <^ n <^ t and t — )• oo 

(4.18) 



Since the term (j2.4p vanishes for /? = 2, for this f3 we can substitute ()4.18p for 5F in (j2.ip 
to include 

log Ef\n;{t,^)) ^^-^ -^ + ^^^^ + -'log^. (4.19) 

i«:n<|t| 

Although there is no exact asymptotic results in this regime, as with (j3.23p the exact 
expansion with (3 = 2 and n fixed is known |3Uj . |23[ eq. (9.89)], and its large t expansion 
is in complete agreement with (j4.19p . 

4.2 V2 in the case n = and general /3 > 

In the case n = 0, P{i){x) is supported on (t, 00) and according to (j4.3p with 61 = 62 = t/2 
substituted in (|2.24p takes on the explicit functional form 

x-t/2 
'^'^^'^ = W^f 

Substituting in (j2.4p and a simple change of variables shows 

T/ 1^ ^3/2 p ^_l/2 X-1/2 

This integral can be evaluated by computer algebra to give 

Adding this to (|2.29p with c = l/vr, then substituting for (5-F in (|2.ip implies 

log Ho; {t, 00)) + (f " ^^'^^^ 

with = {V2 — l/2)/3. However, this contradicts the known [H H] t — >• —00 asymptotic 
form of £^^°f*(0; (t,oo)), which tells us that (ji:2T|) holds with k = V2/3. Thus, as for the 
hard edge, this infinite system implementation of the log-gas formalism breaks down in its 
form of V2- But we observe that the correct value of k follows if we replace V2 by (13.221) . 



19 



4.3 V2 for general 1 <^ n <^ t and general f3 > 

The integration method used by Dyson in [14J to deduce (j3.20p allows us to deduce that 



(4.22) 



We don't give the details since we have already established that this term does not correctly 
reproduce known exact results. 

Instead, we proceed by analogy with the bulk, the corrected hard edge result and the 
corrected n = soft edge result, and thus hypothesize that (|4.22|) should be replaced by 
()3.22p . Doing this, deducing the asymptotic form of V2 using ()4.14p and adding it to ()4.18p . 
then gives 



^ t— 00 

Kn<|t| 



'24 + 3 ' ' 



/3n + 



/3 



+ 



/5 2^ 
2^ + 



1 n 



logltl 



-3/4 



(4.23) 



As for Ej^ {n; (0,t)), in the case n = the form of the term 0(log|t|) in the t — )• —00 
asymptotic expansion of log E^'^^'^{n; (0, t)) is known, and this added to (j4.23p suggests that 
uniformly in n, for n <ti t, the expansion (jl.Sp as recorded in the Introduction holds true. 

Again as for the hard edge case, there is an asymptotic functional relation which (jl.Sp 
must satisfy [22]. This states 



logEf\n;{spt,oo)) ~ 



log Elf^ {pn/2 + p/2 - 1; {s^/^t, 00)) , (4.24) 



where S/3,S4/^ are length scales related by (/3/2)^/'^S/3 = 54/^, and this is indeed exhibited 
by (jl.Sp . It is similarly true that (jl.Sp is consistent with the soft edge analogue of (|3.24p 

m 



(n;(t,oo)) ~ Sr*(n;(t,oo))Sr*(n + l;(t,oo)). 

t— >— 00 



5 Concluding remarks 

An infinite log-gas formalism, due in the bulk to Dyson [14j and independently Fogler and 
Shklovskii [18] , has been applied to the computation of the conditioned soft and hard edge 
gap probabilities E^^^^{n; (t, 00)) and E^^^'^{n; (0, t)). For this purpose the hard and soft 
edges are characterized by their asymptotic densities (jl.6p and (jl.5p . which are taken to be 
the exact profiles of the background densities in the log-gas. The hypothesis (j2.ip asserts 
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that the conditioned gap probabihty is given in terms of the change in free energy resulting 
from the conditioning. In the bulk, it was found in [T3] and [T5] that the first three orders 
in the large t expansion of E^^^^{n] (0, i)), with <S n <S t could be obtained by asserting 
that the change in free energy is due to an electrostatic and entropy term. The latter 
happens to be simply related to the potential drop in going from the bulk of the Coulomb 
fluid, to the region containing the conditioned n eigenvalues (recall (I3.22P above). 

At the hard and soft edges, it is found that following the same ansatz as used in [T3] 
and [18] leads to an inconsistency with known asymptotic results for the n = case of 
^soft^j^. (^^qq)) and E^^'^'^{n; (0, t)). This inconsistency is due to the entropy term. On the 
other hand, if we use instead (j3.22p . which in the bulk is identical with the entropy term, 
consistency with the known n = results is obtained. Making this replacement for general 
<C n <C t allows us to make the asymptotic predictions ()1.7p and (jl.Sp . 

The infinite log-gas approach should be contrasted to Dyson's pioneering strategy [12J, 
applied to the computation of the asymptotic formula (jl.ip for E^^^^{n; (0,t)) with n = 0. 
The calculation consists of two steps: first use the log-gas hypothesis (12. ip to compute the 
large N form of the probability that there are no eigenvalues in the angular interval [0, a] 
in the circular /3-ensemble. This is a large deviation quantity as there would ordinarily be 
on average aN/2'K eigenvalues in this interval. Second scale the variable a in this large 
deviation formula by writing a = 2nt/N to deduce (II. Ij) . It was over 30 years later that 
the log-gas argument was extended to the infinite system [21 [18] , and applied to the study 
of E^"''^(n;(0,t)). 

Large deviation log-gas studies of the soft edge of the Gaussian /3-ensemble were ini- 
tiated by Dean and Majumdar [TJ |8], and by Vivo, Majumdar and Bohigas [35] at the 
soft edge of the Laguerre ensemble. By an appropriate scaling the first term in (I4.2ip 
was deduced. In these works, only the electrostatic contribution Vi was considered. More 
recently [4], using a formalism based on topological recursions pioneered by Eynard and 
collaborators, higher order correction terms to the soft edge gap probability of the Gaussian 
/3-ensemble were obtained (see pTj eq. (20)] for its explicit functional form). Appropriate 
scaling then gave the expansion p.8p in the case n = 0, together with the explicit form of 
the next term in the asymptotic expansion which is a /3-dependent constant (in relation 
to the latter, see also [HlH]). The significance of this from the log-gas perspective is that 
the first two terms of the topological recursion are precisely Vi and V2 from the log-gas 
formalism. In particular, unlike our finding for the infinite log-gas approach at the soft 
edge, the large deviation log-gas computation does correctly give the first two terms in the 
asymptotic expansion. 

This same conclusion holds true at the hard edge. The work of Chen and Manning [5], 
applying what can be viewed as a large deviations log-gas argument at the hard edge of 
the Laguerre ensemble, then performing a double scaling limit to get the asymptotics of 
the hard edge gap probability correctly gives terms the first two (and the third also) terms. 

The question then remains as to why the entropy term V2 in the double scaling limit of 
the large deviation approach leads to correct results for the asymptotic forms of El°^^{n; {t, 00)) 
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and E'^'^'' (n; (0,t)) with n = 0, whereas in the infinite log-gas approach this term must be 
replaced by (|3.22|) . And why in the bulk (j3.22p is an identity for V2 so no modification is 
required. One possible explanation for the breakdown of the form (j2.4p for V2 at the edge 
but not in the bulk relates to the approach to the edge and bulk states from a finite system. 
The log-gas on a circle has uniform density for all N and so approximates the bulk without 
any boundary term corrections. In contrast, the hard and soft edge states for finite A^, as 
distinguished by the density profile, are only local regions of the spectrum. But the entropy 
term V2 (12. 4p is a global quantity, and carries information about the global spectrum which 
is lost if starting from a forever extending edge state as in the infinite log-gas formalism. 
If we accept this, the remaining question, on which we have no immediate answer, is to 
explain the explicit form of the corrected formula for V2, (j3.22p . 
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